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Abstract 

We investigate the time evolution of Hagedorn wavepackets by non-Hermitian quadratic 
Hamiltonians. We state a direct connection between coherent states and Lagrangian 
frames. For the time evolution a multivariate polynomial recursion is derived that de¬ 
scribes the activation of lower lying excited states, a phenomenon unprecedented for Her- 
mitian propagation. Finally we apply the propagation of excited states to the Davies- 
Swanson oscillator. 
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1 Introduction 

In the last two decades considerable interest in non-selfadjoint operators has developed, and 
even the simplest examples have provided phenomena that greatly differ from what is es¬ 
tablished in the familiar Hermitian context. This pronounced deviation is at the core of the 
theory of quantum physical resonances jMoillj and has strongly motivated the research on 
the pseudospectrum of non-selfadjoint operators [TE05j . Our investigation of non-Hermiticity 
will concentrate on the initial value problem 

i edtU{t) = Op [H t ]U(t) , U(0) = Id , 

where e > 0 is a fixed positive parameter and Op [Ht\ is the Weyl quantised operator of the 
quadratic function 

H t (z) = \z ■ H t z , zEIR 2 ”, 

associated with a possibly time-dependent complex symmetric matrix Hf E C 2nx2n . This 
seemingly simple model problem already encorporates several non-Hermitian challenges and 
clearly hints at the behaviour of more general systems in the semiclassical limit e —> 0. 

So far, non-Hermitian harmonic systems have been mostly analysed from the spectral 
point of view or in the context of VT symmetry, see for example mm §3],in av99b| or 
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more recently {GGHS121 IKSTV15] , It has been proven that the condition number of the 
eigenvalues of non-Hermitian harmonic systems grows rapidly with respect to their size |DK04l 
IHenl4] , while spectral asymptotics have been obtained for skew-symmetric perturbations of 
harmonic oscillators as well as for non-selfadjoint system with double characteristics |GGN09[ 
1HPS13) . A complementary line of research IGSI 1. G.SI2. has emphasised the new, unexpected 
geometrical structures that emerge when propagating Gaussian coherent states. Our aim here 
is to extend these geometrical findings to the larger class of Hagedorn wavepackets and to 
add the explicit description of additional non-Hermitian signatures of the dynamics. 

Over decades, the Hermitian time evolution of Gaussian and Hagedorn wavepackets has 
evolved into a very versatile tool with wide-ranging application in many areas. It was realised 
by Hepp and Heller |Hep74[ lHe!75, Hel76 ] that in order to compute the propagation in the 
semiclassical limit one only needs one classical trajectory t e->■ zt through the centre of the 
wavepacket and the linearisation of the classical flow around it, 

St = nD 2 Ht(z t )St , St € c 2nx2n , n = _I 0 d 

where D 2 Ti t (zt) denotes the Hessian matrix of Tit evaluated in zt- This method was widely 
used in applications, in particular in chemistry |Lit86l lYUOOj . More recently Gaussian 
wavepackets have been more systematically used as a tool for the numerical analysis of highly 
oscillatory initial value problems, either within the wide framework of Gaussian beams meth¬ 
ods, see for example LRT13 , or for Hermitian quantum dynamics with Hagedorn wavepackets 
[Lub08i 1GH14] . Let us give a brief overview of the concepts we develop in this work. 

A Gaussian coherent state is parametrised by a phase space point zq G R 2n and a nor¬ 
malised Lagrangian frame, that is a rectangular matrix Zq G C 2nxn satisfying the conditions 

ZqQZq = 0 and Z^QZq = 2ild n . (1.1) 

Assuming zo = 0 for the ease of notation, we introduce the associated lowering operator 

Mz,,) - , i - (-‘f-) , 

and the raising operator A^(Zq) = — A(Zq) as its formal adjoint. The wavepacket cpo(Zo) is 
defined as an element in the kernel of the lowering operator A(Zo), i.e., by A{Zq)^q{Zq) = 0, 
and it is normalised according to 

l|^o(^o)|| 2 = [ \<Po(Z 0 )(x)\ 2 dx = 1. 

JR n 

f This determines <po(Zq) up to a phase factor. The time evolution of the raising and lowering 
operator is easily computed as 

U(t)A(Z 0 )U~\t ) = A(S t Z 0 ) , U(t)A\Z 0 )U-\t) = A\S t Z 0 ) , 

where St denotes the complex conjugate of the complex symplectic matrix St G Sp(n, C). For 
Hermitian dynamics, St is real and StZo = St.Zo is a normalised Lagrangian frame. In the 
non-Hermitian case, however, SfZ q and SfZo both violate the second, normalising, condition 
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of (1.1), and we have to use a Hermitian matrix Nt E C nxn to construct a normalised 


Lagrangian frame Zt := Sf ZoNf with the same range as SfZ q. We then obtain 


U(t)<p 0 (Z 0 ) = t&MZt) 


with ||<£>o(^)|| = 1, where the real-valued gain or loss parameter 


Pt 


tr(G T 1 Irn H r )dr 


is determined by the symplectic metric Gt = e{ZtZ^)VL E Sp(n, R) associated with the 
normalised Lagrangian frame Zt, see also [GS121 §3]. 

The Gaussian wavepacket <po(Zo) is the zeroth element of an orthonormal basis of L 2 (R n ) 
constructed by the repeated application of the components of the raising operator to the 
coherent state. The definition 


<p a (Z 0 ) = A\Z) a cp 0 (Z 0 ) , a E INq , 


is due to Hagedorn |Hag98| , who called the basis elements semiclassical wavepackets, simplify¬ 
ing his earlier construction that was based on a less transparent polynomial recursion |Ha g85| . 
The generalised coherent states { CR121 §4.1], that are obtained by applying a unitary squeez¬ 
ing transformation to the n-fold product of univariate harmonic oscillator eigenfunctions, only 
differ by a phase factor from the semiclassical wavepackets. For our study of non-Hermitian 
dynamics we follow Hagedorn’s ladder approach and use the time evolution of both the raising 
and lowering operators to explicitly describe the propagation of the basis functions. 

Our main new result is the expansion of the time-evolved wavepackets with respect to 
the orthonormal basis cp a (Zt), a E INg, that is parametrised by the normalised Lagrangian 
frame Zf . In the Hermitian situation, we simply have 


U(t)<p a (Z 0 ) = tp a (Zt) , QEljf ■ 


For non-Hermitian dynamics, however, the propagated excited states are utterly different. 
The two Lagrangian frames StZo and StZ$ do not only lose normalisation but also have 
different ranges. Therefore, the dynamics also activate lower lying excited states and we 
obtain 

U(t)<p a (Z 0 ) = e* ^2 a k {t)ip k (Z t ) 

|fc|<H 

with expansion coefficients a,k(t) E C for |fc| < |a|. These coefficients can be explicitly inferred 
from Theorem 4.5, that proves 

u(t)(p a (Zo) = ^=q a (N t A\Z t ))MZ t ) , 


where the multivariate polynomials q a satisfy the recursion relation 


qo(x) = l, q a+ej (x) = Xjq a (x) - ej ■ M t \7q a (x) , j = l,...,n. 
The complex symmetric matrix 

M t = l(S t Z 0 ) T G t (S t Z 0 ) E € nxn 
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governing the recursion is determined by the symplectic metric Gt and the complex flow St- 
In general, the complex symmetric Mf does not have a specific sparsity pattern so that all 
the n dimensions are coupled within the polynomial recursion. 

We have organised the paper as follows: Section [2] develops the symplectic linear al¬ 
gebra of complex Lagrangian subspaces required for the parametrisation of the Hagedorn 
wavepackets. Section [3] constructs coherent states, ladder operators and Hagedorn wavepack- 
ets parametrised by positive Lagrangian frames. Section [4] is the core of our manuscript. 
It analyses the non-Hermitian time evolution of Hagedorn wavepackets, and in particular 
proves our main result Theorem 4.5 Section [5] illustrates our results for the one-dimensional 
Davies-Swanson oscillator. The three appendices summarise elementary facts on Weyl cal¬ 
culus, present a proof of the Riccati equation for the symplectic metric Gt, and discuss basic 
properties of the multivariate polynomials q a , a £ INq. 


2 Lagrangian subspaces 


We start by discussing some symplectic linear algebra with a focus on complex vector spaces 
and complex matrices. We endow the real vector space R 2n with the standard symplectic 
form R 2 " x R 2 ” —> R, (x, y) i —> x ■ liy, using the invertible skew-symmetric matrix 


n = 



-Id„\ 
0 ) 


e R 


2nx2n 


Matrices S £ R 2nx2 n respecting the standard symplectic structure satisfy S T US = H and 
consequently S' -1 = Q T S T Ll. They are called symplectic and constitute the symplectic group 
Sp(n, R), see also [MB'98 : §1.2]. Writing a symplectic matrix as S = (U, V) with U,V € R 2nxri , 
the complex rectangular matrix Z = U — iV £ C 2nxn satisfies 

Z T nZ = 0, Z*VlZ = 2ild n , (2.1) 

—T 

where Z* = Z denotes the Hermitian adjoint. We see from the first property of Z that all 
vectors l, l' £ range Z satisfy 

i-m' = o. 

Such vectors are called skew-orthogonal, and a subspace L C C n © C n is called isotropic, if 
all vectors in L are skew-orthogonal to each other. Moreover, L is called Lagrangian, if it 
is isotropic and has dimension n, which is the maximal dimension an isotropic subspace can 
have (by the non-degeneracy of H). From the second property of the matrix Z, we see that 
all vectors l £ range Z \ |0| satisfy 

i(nz) -l > 0 . 

That is, the quadratic form 

h{z, z') := i(fif) • z' = \z ■ Vi T z!, z, z' £ C n © C n , 
is positive on range Z. Such a Lagrangian subspace is called positive. 

Remark 2.1. In the literature, ILubOStf . the choice Z = U + iV is more common, but for our 
purposes the complex conjugate is more natural, since it matches with Hagedorn’s notation, 
\Hag9$j. However, the main results hold true for both definitions. 
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2.1 Lagrangian frames 


Rectangular matrices satisfying conditions (2.1) are convenient tools when working with La¬ 
grangian subspaces. In particular, the normalisation condition Z*OZ = 2ild n will be crucial 
later on when studying the effects of non-Hermitian dynamics. 


Definition 2.2 (Lagrangian frame). We say that a matrix Z £ C 2nxn i s isotropic, if 


z T oz = o, 


and, it is called normalised, if 

Z*OZ = 2ild n . 

An isotropic matrix of rank n is called a Lagrangian frame. 

As indicated before, normalised Lagrangian frames are in one-to-one correspondence with 
symplectic matrices: Writing S £ Sp(n, R) as S = (U, V ) with U, V £ R 2nxn , then Z = U—iV 
is isotropic and normalised. Vice versa, if Z £ C 2nxn is a normalised Lagrangian frame, then 
S = (Re(Z), — Im (Z)) is symplectic. For a positive Lagrangian subspace L C C n © C n there 
are plenty of normalised Lagrangian frames spanning L. Denoting by 

F n (L) = {Z £ C 2nxn ; range Z = L, Z*OZ = 2ild n } 

the set of normalised Lagrangian frames spanning L, we observe that all its elements are 
related by unitary transformations. Indeed, since any Zq,Z\ £ F n (L ) have the same range, 
there exists an invertible matrix C £ C” xn so that Z\ = ZqC, and the normalisation requires 
that C is unitary. 

2.2 Orthogonal projections 

The complex conjugate L of a positive Lagrangian subspace L C C n © C n is Lagrangian, too, 
and all vectors l £ L \ {0} satisfy 

h{i,i) = ±I-n T i < o , 

so that L is called a negative Lagrangian. It is clear that L n L = {0}, because if l £ L n L, 
then l is real, and hence h(l,l) = i l ■ Q.l/2 = 0, so that l = 0. Therefore, 

<D n © C n = L® l . 

This decomposition of C n © C n is orthogonal in the sense that 

h{l. I') = 0 for all l £ L, l' £ L . 

Proposition 2.3 (Projections). Let L C C n © C n he a positive Lagrangian and Z £ F n (L). 
Then, 

7 t l = \ZZ*O t and ir L = -±ZZ T 0 T 
are the orthogonal projections onto L and L, respectively, that is, 

(i) | l= Id 2 n, 1 1= 0 and ir L \ L = Id 2n , \ L = 0 , 

(ii) 7r 2 = 7 tl and 7r| = 7r^ , 
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(Hi) h(-KLZ , z') = h(z, ttlz') and hQTpz , z') = h(z, 'Kpz') for all z, z' £ C n © C n , 

Proof. To prove ttl\l= Ichn and ttl \ f= 0, we observe 

7 x L z = \zz*vl t z = z , 7 t l ~z = ±zz*n T z = o . 

The other properties of 717 , and 7 xp are also proved by short calculations using that Z is 
isotropic and normalised. □ 

2.3 Siegel half space 

A large set of Lagrangian subspaces can be naturally parametrised by complex symmetric 
matrices. If the Lagrangian is positive or negative, then we encounter complex symmetric 
matrices with positive or negative definite imaginary part, that is, elements of the upper or 
lower Siegel half space. 

Lemma 2.4 (Siegel half space). Assume that L C C n © C n is a Lagrangian subspace so that 
the projection C n © C n —> C n , (p, q) i->- p is non-singular on L. Then there exists a unique 
symmetric B G C nxn such that 


L = {(Bq,q)- q eC n } . 

The matrix B can be written as B = PQ -1 , where P,Q G C nxn are the components of any 
Lagrangian frame Z G C 2nxn spanning L, that is, 


Z = 



range Z = L . 


Furthermore, L is positive (negative) if and only iflmB is positive (negative) definite. 

Proof. That the projection of L to C n is non-singular means that there is a function f(q) 
such that L = {{f(q),q ); q G C n } and since L is linear / has to be of the form f(q) = Bq 
for a uniquely determined matrix B G <D nxn . Now let us denote Ib(q) = ( Bq,q ) for q G C n . 
Since L is isotropic, we must have 

0 = l B (q) ■ m B (q') = q ■ (B - B T )q' 

for all q, q’ G C n , hence B = B 1 . If Z = ( P,Q ) and Z\ = (P\]Qi) are Lagrangian frames 
spanning L , then there is an invertible matrix C G C nxn with Z\ = ZC , so that 

PiQf 1 = PQ _1 = B . 


Furthermore, 

h(l B {q),lB(q)) = ^lB{q)*^ T lB(q) = h,q*{B - B)q = q*lm.Bq 
for all q G C n , so that L is positive (negative) if and only if ImB is positive (negative). □ 


6 



2.4 Metric and complex structure 

The Hermitian squares of normalised Lagrangian frames have been useful for writing projec¬ 
tions on Lagrangian subspaces. We now examine their real and imaginary parts to see more 
of their geometric information unfolding. 

Proposition 2.5 (Hermitian square). Let Z £ C 2n,xn be a normalised Lagrangian frame. 
Then, 

ZZ* = R e{ZZ*)-in , 

where R e(ZZ*) £ Sp(n, R) is a real symmetric, positive definite, symplectic matrix. In par¬ 
ticular, R e(ZZ*)~ 1 = fl T Re(ZZ*)Ll.. Moreover, 

R e(ZZ*)QZ = i Z , Re{ZZ*)LlZ = -i Z , 

so that (R e(ZZ*)fl) 2 = —M 2 n . 

Proof. Writing ttl + irp = Id. 2 n in terms of Z, we obtain — Im (ZZ*)Ll T = Id 2 n- Hence, 
Im (ZZ*) = Q T . This implies symplecticity of the real part, since 

Re(ZZ*) T nRe(ZZ*) = \{ZZ T + ZZ*)Vl(ZZ T + ZZ*) 

= \(2iZZ* - 2iZZ T ) = -Im(ZZ*) = Q . 

Checking positive definiteness, we see 

2 ■ Re(ZZ*)z = \z ■ {ZZ*z + ZZ T z) = \Z*z\ > 0 

for all z £ H 2n . If Z*z = 0, then ZZ*z = 0 and Im(ZZ*)z = 0, which means z = 0. Finally 
we compute iR e(ZZ*')QZ = |( ZZ* + ZZ T )LIZ = —Z. □ 

We have already observed that two normalised Lagrangian frames Zq,Z\ £ F n (L) are 
related by a unitary matrix C £ C” xn with Z\ = ZqC. Therefore the Hermitian squares 
ZqZq = Z\Z\ are the same and can be used for defining two key signatures of the La¬ 
grangian L. 

Definition 2.6 (Metric & complex structure). Let L C t"0 C n be a positive Lagrangian 
subspace and Z £ F n (L). 

(i) We call the symmetric, positive definite, symplectic matrix 

G = n T R e(ZZ*)n 


the symplectic metric of L. 

(ii) We call the symplectic matrix 

J = -QG 

with J 2 = —Id 2 n the complex structure of L. 

The complex structure J is a symplectic matrix so that LlJ is symmetric and positive def¬ 
inite. Such complex structures are called H-compatible. That positive Lagrangian subspaces 
and H-compatible complex structures are isomorphic to each other, has been observed and 
proven in |GS12l Lemma 2.3]. The complex structure can also be used for concisely writing 
the orthogonal projections. 
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Corollary 2.7 (Orthogonal projections). Let L C C n © C n be a positive Lagrangian and 
J E Sp(n, R) its complex structure. Then the orthogonal projections on L and L can be 
written as 

tt l = |(Id 2n + iJ) , 7 t l = \(lchn - i J) ■ 


Proof. Propositions T3 and |2.5| yield 

ir L = i ZZ*n T = i Re(ZZ*)U T + ild 2n = |(Id 2n + iJ) 


□ 


We can construct a normalised Lagrangian frame from the eigenvectors of the matrix 
G E Sp(n, R) representing the symplectic metric. To this end recall the basic structure of the 
spectral decomposition of a positive definite symplectic matrix G. 

Lemma 2.8 (Spectrum of symplectic metric). Suppose G E Sp(n,R) is symmetric and 
positive, then there exists a basis u \,..., u n , v ±,..., v n E R 2ri such that 

Gu]$ — ^k'U'k ? Gvk — A fc Vk , 


where Xk > 1 , for k = 1 , ■ ■ • ,n, and for all j, k = 1 , • • ■ n we have uyUk = Vj-Vk = Vj-Lhik = 5jk 
and Uj ■ Lluk = Vj ■ Llvk = 0 . 

Proof. This result is in principal well known, see e.g., [MS981 Lemma 2.42] for a similar 
statement, but it is hard to locate this exact form of it, so let us indicate the basic idea. 
Since G is symplectic we have GLl = LlG^ 1 , and since G is symmetric there exists a basis of 
eigenvectors. Now let u\ be an eigenvector with eigenvalue Ai > 0, then v\ := Vlu\ satisfies 
Gv i = GLlu\ = LlG~ 1 u i = A^Om = A]” 1 ^, and hence is an eigenvector with eigenvalue 
Ajj 1 . So we can assume Ai > 1 , and u\ ■ v\ = 0 since G is symmetric, and if we normalise 
u\ as ui ■ u\ = 1 then v\ ■ Qui = 1. Let Vi be the span of u\,v\, then V f 1 = Vj n , where 
Vj := {u E V : v T Qu = 0 Vu E Vi}, this follows since with v\ = Jlwi and u\ = —Aim the 
conditions v ■ u\ = 0 and v ■ v\ = are equivalent to v ■ Qui = 0 and v ■ Llu\ = 0. Therefore 
R 2 n _ \/j 0 an( t yf- j s symplectic and invariant under G, hence we can repeat the previous 
step in Vi~ and arrive after k steps at a basis with the properties claimed in the lemma. □ 


Lemma 2.9 (Normalised Lagrangian frame). Let G E Sp(n, R) be symmetric and positive 
definite. Consider an eigenbasis u \,..., u n , v \,..., v n E R 2n ofG as described above in Lemma 
13.81 and denote 

Ik ■— i\/AfcUfc , k — 1 ,...,n . 

Then, the matrix Z E C 2nxn with column vectors l \,..., l n is a normalised Lagrangian frame 
so that G = fl r R e(ZZ*)Ll. 


, v n from Lemma 2.8 we find lj-Lllk = 


Proof. Using the properties of the basis u\, ■ ■ ■ u n , v\, 

0 and lj-Lllk = 2i 5jk, so Z is a normalised Lagrangian frame. Furthermore, again using Lemma 
2.81 we obtain 


Re(ZZ*) = J2Mhll) = £ (& k uZ + A kVkvl) = G~ l . 

k=1 k=1 


□ 
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3 Raising and lowering operators 


Coherent states can be characterised by their lowering operators, or annihilators. These are 
operators with linear symbols, so let us briefly define them and review some of their properties. 
We will denote 



where (fnf)(x) = —ieV x tjj(x) is the momentum operator and (qif)(x) = xip(x) the position 
operator. 

Definition 3.1 (Ladder operators). Let I G t n $ C n , then we will set 

A(l):=^=l-Llz , A\l) := -A(J) , (3.1) 

V 2e 

A(l) is called a lowering operator, while A\l) is called a raising operator. 

The following properties are important but easy to prove. 

Lemma 3.2 (Commutator relations). We have for all l, l' E C n © C n 

(i) [A(l),A(l')\ = -±l-Lll' , 

(ii) [A(l),A\l')] = \l-LlV = h(l',l) . 

A*(l) is (formally) the adjoint operator of A(l). 

Proof. We use the phase space gradient V = V Pi9 . Basic Weyl calculus, see appendix [Aj 
implies for any symbol b 

[A(l), Op[6]] = ieOp[VA(Z) • LlVb] = sJJJ 7 ! Op[Z • V6] , 

since 

VA(l) • Ll\7b = -^LlU • fiVft = =1 • V6 . 

\[2e V2e 

Choosing b = ^=Z / • Liz, we have V6 = — -^=Lll' and this gives us (z). Part (ii) follows with 
choosing b = — -AJf ■ Liz instead. □ 

We see in particular from the first property that we can create a set of commuting lowering 
operators if we choose a set of Vs which are skew-orthogonal to each other. Moreover, a 
Lagrangian subspace parametrises a maximal family of commuting lowering A(lf s. Following 
Hagedorn |Hag98| , we combine them as an operator vector. 

Definition 3.3 (Ladder vectors). For an isotropic matrix Z 6 C 2nxn with columns l \,..., l n 
we will denote by A(Z) and A^(Z), the vectors of annihilation and creation operators, respec¬ 
tively, 


A(Z) := (A(h),-.. ,A(l n )f = -^=Z T Llz , 
A f (Z) := (A\h),... ,A\l n )f = ^=Z*Qz . 
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For any multi-index a G INq , we set 


A a (Z) := A(l 1 ) a 'A(l 2 ) a *---A(l n ) a ” , 

Al(Z) := A\hrA\l 2 r-..A\l 71 ) an ■ 

Since all the columns of an isotropic matrix are mutually skew-orthogonal, all the com¬ 
ponents of the annihilation vector A(Z) commute. The same is true for the creation vec¬ 
tor A^(Z). Therefore, the operator products A a (Z) and Af(Z) do not depend on the ordering 
of their individual factors. 

Remark 3.4 (Hagedorn’s parametrisation). The ladder parametrisation coincides with the 
original one of Hagedorn \Hag98\ [ /. Considering matrices A,B€ C nxn with A T B — B r A = 0 
and A*B + B* A = 2Id n , he sets 

An^ g (A,B) = ~^=(B T q + iA T p) , 

<(AB| = p) . 

We can write 

Au ag (A, B) = -J= (-i B T q + A T p) = -j=Z T nz with Z = f 1 ^) , 
and quickly convince ourselves that Z is isotropic and normalised as well. 

3.1 Coherent states 

Coherent states emerge as a joint eigenfunction with eigenvalue 0 of a family of commuting 
operators parametrised by a Lagrangian subspace L C C n © C n . We set 

I(L) := {(p G V'(R n ) ; A(l)tp = 0 for all l E L} , 


and observe that 

ip G I{L ) if and only if A(Z)<p = 0 

for any Lagrangian frame Z G F n (L). The following characterisation of I(L) is quite standard. 
For instance one can find a similar statement in [Hor951 Proposition 5.1], but let us sketch 
the proof to elucidate how the Lagrangian property implies the Gaussian form. 

Proposition 3.5. Consider a Lagrangian subspace L = {(Bq,q ); q G C n } parametrised by a 
symmetric matrix B G C nx?l . Then, every element in I(L) is of the form 


p(x)=ce^ x - Bx , 

for some constant c G C. Furthermore, L is positive if and only if I(L) C L 2 (R n ). 


Proof. As in the proof of Lemma 2.4 we denote Ib(x) = (Bx, x) for x G C n . Let l G L. 


A(l)ee xBx / 2 = —l ■ ni B (x)e^ 

v 2 £ 




= 0 


(3.2) 

Then, 
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using that l B {x) £ L implies l ■ Lll B (x) = 0. Hence e^ x ' Bx £ I(L), and ei x ' Bx £ L 2 (H n ) if 
and only if Im5 > 0, which is equivalent to the positivity of L. To show uniqueness we use 


that 


V2e 


e^ xBx fa e~± x ' Bx 


A(l B (ei)) . 


If (p £ /(L), then we find 


di(e-&- Bx tp(*)) = e-^ xBx J^A(l B ( ei )Mx) = 0 

for i = 1,... ,n, therefore <p(x) = ceFi x ' Bx for some c £ C. Here we used that L has dimension 
n. □ 

Hagedorn’s raising and lowering operators |Hag98| originate from his earlier parametrisa- 
tion of coherent states |Hag85| , which can be conveniently expressed in terms of Lagrangian 
frames. 


Lemma 3.6 (Coherent states). Let L C C” © C n be a positive Lagrangian and consider a 
Lagrangian frame Z £ C 2nxri spanning L. Define P,Q£ C nxn by 



Then, Q and P are invertible and 

<p 0 (Z;x) := {TT£)- n/i (detQ)-^ 2 e^ x - p Q~ lx € I(L) . (3.3) 

Furthermore, Z is a normalised Lagrangian frame if and only if 

\\To(Z)\\ 2 = [ \ipo(Z] x)\ 2 dx = 1. 

J E" 

If C £ C nxn is non-degenerate then 

MZC) = (det cy^Mz) • ( 3 . 4 ) 


Proof. Rewriting positivity of the Lagrangian L in terms of P and Q gives 

±z*nz = \{P*Q-Q*P ) >0 . 

Hence, ^((Py)*(Qy) — ( Qy)*{Py )) > o for all y £ C n so that P and Q are invertible. Then 
Lemma 2.4 and Proposition 3.5 imply that the Gaussian wave packet of (3.3) is an element 
of I(L). The normalisation of Z is equivalent to ^(P*Q — Q*P ) = I, and multiplying from 
the left by Q * _1 and from the right with Q -1 gives 

±(PQ- 1 -Q*~ 1 P*) = (QQ 


*\-l 


which is the same as 


ImiPQ- 1 ) = (QQ*) 
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This implies that po(Z) is normalised, since then 

J\(po(Z; x)\ 2 dx = (7re) _n,/2 |det Q| _1 J e~£ x ' lmP ® ±x dx 


= |detQ| 1 (detImPQ 1 ) 


-n-i /2 _ 


= 1 . 


The relation between the states < po(Zi ) with Z\ = ZC and po(Z) follows by observing that 
Pi = PC and Qi = QC , hence QiP^ 1 = QP~ X and detQi = detQdetC. □ 


Notice that (3.3) defines <po(Z;x) only up to a phase, because we have not specified the 
branch of the square root of det Q. In practice it will typically be determined by continuity 
requirements. 

3.2 Orthonormal basis sets 

Let L C C n © C n be Lagrangian. Applying the operators (l) multiple times to an element 
in I(L ) will be used to create a basis. To see the basic idea assume L is positive and <po £ I(L) 
has norm one, 

12 = (ipo,Po) = / Po(x)vo (x)dx = 1 . 


|<Po | 


|R" 


Then we can use the relation 

A{l)A\t) = [A{1 ), A\l')] + A\l)A(l') = h(l', l ) + A\l)A{l') 

to obtain 

(A\l)w,A\l')w) = {po, A(l)A\l')<p 0 ) 

= h(l',l)(p 0 ,Po) + (cpo,A\l)A(l')cp 0 ) 

= h{l',l ) , 

where we have used that A(l')ipo = 0. So if h(l',l ) = 0, then the states A^(l)ipo and A\1')lpq 
will be orthogonal to each other. It is easy to check that they are both orthogonal to (/?o and 
that ||At(0 po|| = 1 if h(l,l) = 1. Iterating this construction yields an orthonormal basis. 

Theorem 3.7 (Orthonormal basis). Let L C C n © C n be a positive Lagrangian subspace and 
Z £ P n (L). Then for any normalized ipo £ I(L) the set 


<Pa(Z) := —=Al(Z)(p 0 , 


a £ IN" 


o > 


(3.5) 


is an orthonormal basis o/L 2 (H" 


This result is due to Hagedorn [Hag98, Theorem 3.3]: The normalisation and orthogonality 
follows from commutator arguments similar to the simple case we discussed. Completeness 
can be derived from the fact that the functions p a {Z ) are the eigenfunctions of the number 
operator 

N(Z) = A(Z) ■ A\Z) = A{h)A\h) + • • • + A(l n )A\l n ) , 

which is selfadjoint and has a complete basis of eigenfunctions due to the following Lemma. 
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Lemma 3.8 (Number operator). Let L C C" ffi C” be a positive Lagrangian subspace. Let 
Z E F n (L) and G E Sp(n, H) be the symplectic metric of L. Then we can write N(Z) = 
A(Z ) • A\Z) as Weyl operator N(Z) = Op[z^] with symbol 


n(z) = — (z ■ Gz + ne) , z E R 2n . 

2e 

Proof. By Lemma |A.1[ 

A(Z) ■ A\Z) = A- (jg . Lt T ZZ*Llz + )|(— 2in)) = ^ (2 ■ Gz + e?r) 


where we have also used that Proposition 2.5 implies ZZ*Q = G — iLl 


□ 


Returning to the previous remark, by Lemma 3.8, N(Z) is the Weyl quantisation of a 


positive definite quadratic form. By symplectic classification of quadratic forms, see [Hor94l 
Theorem 21.5.3], such a form is symplectically equivalent to a sum of harmonic oscillators, 
and using the quantisation of linear symplectic transformations as metaplectic operators, see 
}CR,12l §2.1.1], N(Z) is therefore unitary equivalent to a sum of standard harmonic oscillators. 

The orthogonality and normalisation of the basis functions <p a (Z) depend on the normal¬ 
isation of the matrix Z. Let us examine the creation process with parameter matrix ZC , 
where C is non-degenerate. Then, (po(ZC) E /(L). However, the next creation step provides 
orthogonality if and only if C is unitary, since 


(« ■p ej (ZC),<p eh (ZC )) = hfZCej, ZCe k ) = e T k C*Cej . 


Let us expand ipp(CZ ), a member of the possibly non-orthogonal function set, with respect 
to the orthonormal basis cp a (Z), a E INq. 

Theorem 3.9 (Expansion coefficients). Assume Z E C 2nxn is isotropic and normalised, and 
let C E C nxn be non-degenerate, then for all a,f3 E INq 

W ZC),^(Z)) = Tpb £ (3,6) 

(aetoj2 Agm ( a ^ 


where we denote 

m(at,/3) := |a = (Ay) E 1N 0 ; X^i=i ^ij = Pj > Ylj= l Ajj = a «} i 

as well as A! := flij A ij\ and C A := Ylij . 

Proof. We have Af{ZC) = C*A^(Z) by dehnition of A^(Z). So if we write y = A^(ZC) and 
x = A^(Z), then we have y = C*x and 

^ (ZC) = ^ y ^ o(ZC) = ^(detC) 1 / 2 ^ 0 ^ ' 


To proceed we have to expand y ;3 in terms of x. Since y* = c, • x, where c t is the Tth column 
vector of C, we obtain = (ci • x)^ 1 • • • (c n • x)^ 71 , and for the individual terms we use the 
multinomial expansion, 


o* = £ § 


AH 


c}V J 


IA 
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where the A J E IN," are multi-indices and for any vector x = (ari, • • • ,x n ) T E C n we set 


A! 


x A ' := n 4 x ? ; ‘ • Multiplying the terms for different j gives 


y" = 


£ 


t© 1 Fr A"„A 1 +-+A" 

C 1 c n x 


|A J |=/3j,j=l,...,n 


A 1 ! • • • A n ! 


Therefore we found 


< P p(ZC) = (detC)~ 1 / 2 £ 

|A 3 'N3,i=l,-,n 


^rrr~>pu c i '" c « aA ^ 1 h i~ ■^ ti )^a i -i—■+a"(-^) 


and taking the overlap with ip a (Z, z ) and using orthogonality gives 

{<pp{ZC),y a {Z)) = (detC)- 1 / 2 £ ■ ■ ■ cf VriS a , x i + ... +X n . 

|A J '|=ft,i=l,--- ,n 

If we introduce the matrix A with columns A 1 ,..., A n , then this formula can be rewritten as 
in the statement. □ 


3.3 Phase space centers 


Sofar we have focused on positive Lagrangian subspaces L and Lagrangian frames Z E F n (L ) 
and have discussed coherent states ipo (Z) centered at the phase space origin. Now we extend 
this framework to formal complex centers z E C n © C n . This generalisation is motivated by 
the choice of complex Hamiltonians. To give a physically meaningful interpretation of the 
associated position and momentum further investigation is needed. 


Definition 3.10 (Centered ladders). For l,z E C n © C n we define the ladder operators 

A(l, z) := ~j=l ■ H(£ — z) , A\l, z ) := — A(l ., z) . 

y/2s 

We note that A(l , 0) = A(l) and A^(l, 0) = A^(l). 


Adding a constant to an operator does not change its commutation properties so that 
Lemma 3.2 also applies to A(l,z ) with z 0, and each of the ladder vectors 

A(Z, z ) := -L Z T n(z - z) , A\Z, z) := - z) , 

y2 £ v 2e 


has commuting components, if Z E C 2nxn is an isotropic matrix. One can change the center 
of ladder operators by conjugating with the (Heisenberg-Weyl) translation operator 


T(z) =exp(-^-P5) , z = (p,q)€ C n ©C n , 

that acts as (T(z)ijj)(x) = e^ p '^ x ~^ qS, ip(x — q) on square integrable functions E L 2 (IR, n ), 
which have a well-defined extension to C n . Indeed, it follows easily from the definition that 
T(w)zT(w )~ 1 = z — w, which directly yields 

T(w)A(l, z)T(w)~ l = A(l, z + w) (3-7) 


for all w,z E C n © C n . Therefore, all the previous results can be translated away from the 
origin. We have: 
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Theorem 3.11 (Orthonormal basis). Let L C C n © C n be a positive Lagrangian subspace 
and Z E F n {L). Let z = ( p, q) E C n © C n . Then every element in 

I(L , 2) := {</? E P 7 (R”); j4(Z, z)<£> = 0 for all l E L} 

is a constant multiple of the normalised coherent state 

p 0 (Z, z; x) := ( 7 re)- n/4 (det Q)- 1 /2 e ^0- g )-PQ- 1 (^-<?)+ip-0-9) ; (3.8) 


and the set 


tp a (Z,z) := -?=A f a (Z,z)<po(Z,z) , a E Mq , 
v a! 


is an orthonormal basis of L 2 (M n ). 


Proof. We use the translation property (|3.7|) to lift the results of Proposition 3.5 and Theo¬ 
rem 


3.7, that apply for 


I(L, 0) = I(L) and <p a (Z, 0) = <p a (Z) , 


to general centers z / 0. □ 

It turns out that one can always reduce to the case with real center z. To understand why 
let us ask which conditions on z, w E C n © C n must hold so that I(L, z) = I(L, w). In terms 
of the annihilation operators this means that for all l E L and ip E I(L, z) 

A(l, z)tp = A(l, w)<p . 


Since 


A(l, z) — A(l, w) 


V2e 


■ Ll{w 


z) , 


this is equivalent to the condition that l ■ — w) = 0 for all l E L. So z 

orthogonal to L, but since L is Lagrangian this means that 


w has to be skew 


z — w E L . 


So any two complex centers whose difference is in L define the same ladder operators, and 
we just have to find v E C n © C n so that w = z + ttlv is real. Now we use Corollary 2.7 and 
write 7 tlv = \{y + i Jv). Then we immediately see that v = — 2ilmz provides the real center 


w = Re z + J Im z 


w — z E L . 


In summary we have reproduced the linear algebra part of [GS12L Theorem 2.1] that also 
provides the associated coherent states. 

Theorem 3.12 (Real centers). Let L C C n © C n be a positive Lagrangian and Z E F n (L). 
Let J E Sp(n, 1R) be the complex structure of L and define Pj : C n ffi C" —> R n © R ra , 
Pj(z) = Re z + J Im 2 . Then, for any z = (p, g) 6 C" ffi C n 

A(Z, z) = A(Z, Pj(z)) , I(L, z) = I(L, Pj(z )) 


and the coherent states are related by 

To(Z,z) = et^+v)iy-q)^ (z,Pj{z)) , (rj,y) = Pj(z ) . 
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4 Time evolution 


We will now explore how the ladder operators, coherent states, and the associated basis behave 
when we propagate them in time according to a non-Hermitian operator. Let Ht £ (C 2nx2n 
be a symmetric matrix, which depends continuously on t, and denote by Ht = Op [Ht] the 
Weyl quantisation of the quadratic function 

H t (z) = \z ■ H t z , z £ R 2ri . 

We are interested in the Schrodinger equation 

i ed t ip{t) =H t r ip(t) (4.1) 


and the corresponding time evolution operator U ( t ) defined by 

i ed t U(t) = H t U(t) , 17(0) = Id . 

If Ht is a real symmetric matrix, then we are in the standard setting. Ht is a self-adjoint 
operator on some dense domain of L 2 (R n ) and U(t) defines for all t £ R a unitary time 
evolution, see for example jCR,121 §3.1]. In the time-independent case with ImH < 0, operator 
theory allows to define the evolution U(t) as a contraction semigroup on L 2 (R n ), see the 
general Mehler formulas in [Hor95|, in particular the discussion before [Hdr951 Theorem 4.2], 
In the more general case that ImHt depends continuously on time and Im Ht < 0 the well- 
posedness of the problem (4.1) in L 2 (R n ), and hence the existence of U(t), follows then from 


the time-independent case and general results by Kato on hyperbolic evolution systems, see 
[ Paz83l §5.3]. However, if Irn H > 0, then U(t) might cease to be well-defined after some finite 
time T > 0, as shown by the examples in Section [5] or |GS12j . Our aim here is to explore the 
non-unitary evolution of ladder operators, coherent states and excited states on time intervals 
for which U(t) is a well-defined operator on L 2 (R n ). We first investigate how a non-vanishing 
imaginary part changes the geometrical structure. 


4.1 Metriplectic structure 

We decompose the Hamiltonian function Ht.(z) into its real and imaginary part and first 
consider the Schrodinger equation for the real part, 


\edt'ib{t) = Op[Re Ht]ip(t) • 


(4.2) 


Since Op[Re7^] is a self-adjoint operator, the conventional Schrodinger equation (4.2) can be 
reformulated as the Hamiltonian equation 


dt'ipit) = X £t (ip(t)) , 


where X £t denotes the Hamiltonian vector field for the energy function 

£tW ■= (ip,iKij;) , € H 2 (R n ) , 

with K = ^ Op[R,e7/t], see for example [MR991 Corollary 2.5.2], Indeed, one equips the 
complex Hilbert space L 2 (R, n ) with the symplectic form 

■= 2Im((p,i/}) , ip,ip £ L 2 (R n ) , 
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and computes for the derivative of the energy 


= i™s(£tW , + M- f t(V0) 

= (- 0 , + (<£>, iKip) 

= —i(K'if>,ip) + i(K'il>,ip)=uj(K'i/j,ip) 

for all ip, £ H 2 { R n ), so that 

E Op[ReHt]^ = X £t {^) . 

Now let us consider the more general case of non-Hermitian time evolution, which is not 
captured by symplecticity alone but requires additional metric structure. We set 

g(<P,4>) : = 2Re{(p,tp) , p, ip £ L 2 (R n ) , (4.3) 

and observe that g is a symmetric and positive definite IR-bilinear form. This metric defines 
the gradient flow contribution generated by the imaginary part. Indeed, setting 

W) := (i>, l Op[Im n t ]ip) , € H 2 (R n ) , 


we have 


= (Vh l Op[Im 'Htjtp) + (p, l Op[Im = g{\ Op[Im n t ]^, p) 

for all p, ij) £ H 2 (H n ), since Opflm'Rt] is self-adjoint. In summary, we can rewrite the non- 
Hermitian Schrodinger equation (4.1) as 


dfip(t) = X £t (ip(t)) + grad X t {^(t)) 


(4.4) 


and we note that such an additive combination of Hamiltonian and gradient structure defines 
a metriplectic system in the sense of IBMR131 §15.4.1], if additional compatibility conditions 
on the energies £t(ip) and -Tt(V0 are satisfied. 

In the following we will see how a similar metriplectic structure emerges as well in the 
semiclassical limit of the propagation of coherent and excited states. 


4.2 Ladder evolution 

Let St £ C 2nx2n be the matrix defined as the solution to 

S t = SlH t St , S 0 = Id 2n , (4.5) 

for some time-interval [0, T[. It is easy to check that St is a complex symplectic matrix, i.e., 

sf ns t = n , 

and if Ht = H does not depend on time, then St = exp(ffliL) exists for all t £ R. If Ht is a 
real matrix, then St will be a real symplectic matrix. Otherwise, the matrix St is complex. 
We first examine the dynamics of the ladder operators with initial center at the origin. 
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Lemma 4.1 (Ladder evolution). We have for all l E C n © C n 


U(t)A{l)U- 1 (t) = A(S t l) , U(t)A f {l)TJ-\t) = A f (S t l) ■ 


Proof. Let A t := U(t)A(l)U 1 (t), then A t satisfies 

ied t A t = [Hu A t ] . 


Since the symbol at of A t is 


linear, Equation (|A.l ) yields 


dt&t - ifHtyOt} , 


where {/, gj = Vf-fiVg denotes the Poisson bracket. By (4.5) 
equation is given by at(z) = ao(Sf 1 z) and thus, we find 


the solution of this differential 


at(z) 


-$= i ■ nsf l z = -^= s t i ■ Sf T nsf'z = ~^= s t i ■ n z . 

f s/Ye 1 1 VTe 


□ 


The previous Lemma implies for any isotropic matrix Zq E C 2nxn 

U(t)A(Zo)U~ 1 (t) = A(S t Zo ) , U(t)A\Z 0 )U-\t) = A\S t Z 0 ) , 

and we observe that both matrices StZ o and StZo inherit isotropy, since St and St are sym- 
plectic. However, even if Zq is normalised, neither StZo nor St.Zo need to be normalised, since 
in general 

s^s t + n. 

Furthermore, the raising operator is no more the adjoint of the lowering operator, 

A\StZ 0 )^A*{S t Z 0 ) , 

and if Lq = range Zq is the initial Lagrangian subspace, then A(StZo) and A'(StZo) belong to 
the different Lagrangian subspaces StL^ and StLo, respectively. Only if Hf is a real matrix, 
then St.Zo stays normalised, while both the raising and the lowering operators are adjoint to 
each other and belong to the same Lagrangian subspace StLo. 

4.3 Coherent state propagation 

Let us next consider the evolution of coherent states on time intervals [0, T[ so that 

L t := S t L 0 

is a positive Lagrangian subspace. Such intervals exist by continuity of t *->• St, if the initial 
Lagrangian Lq is positive. If we propagate a normalised Lagrangian frame Zq E F n (Lo) by the 
complex flow matrix St, then St.Zo is in general not normalised and the associated coherent 
state tpo(StZo) is not normalised either. Hence, we look for a normalised replacement of St.Zo- 
Since Lt is positive by assumption, the matrix 

N t ■.= (^( S t ZoTSl(S t Zo ))“ V2 
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is well-defined, in particular Hermitian and positive definite, so that 

Z t ■= S t Z 0 N t G F n (L t ) . 

We note that for every normalised Lagrangian frame Wt G F n (Lt) there exists a complex, 
invertible matrix C\ G C nx ” such that Wt = StZoCt- The Lagrangian frame Zt is a particular 
one in the sense that Nt is Hermitian and positive definite. This property will allow us to 
explicity write the time evolved coherent state U(t)ipo(Zo) in terms of a normalised coherent 
state accompanied by a positive loss or gain factor. 

Proposition 4.2 (Coherent state evolution). Let Lq C C n © C n and Lt = be positive 
Lagrangian subspaces for t G [0,T[. Let Gt G Sp(n, R) be the symplectic metric of L t and 
consider Z t G F n (L t ) so that Z t = StZe>N t for a Hermitian positive definite matrix N t G C" xn . 
If the coherent state tpo(Zo) is given by ( |3.3[ ), then 

U(t)<p 0 (Zo) = MStZ 0 ) = eP*<p 0 (Zt) , t G [0, T[ , 

with 

Pt = ~i [ ^(Gf 1 Im H T )dr . 

Jo 

Proof. The proof for real matrices Ht is well known and goes back to Hagedorn [Hag80| . It 
can be extended without any changes to the complex case to obtain 


U(t)<po(Z Q ) = MStZo) , 

see the proof of Proposition |4.8| later on. Switching to the normalised Lagrangian frame Zt, 
Lemma 3.6 implies 

U(t)ip 0 (Z 0 ) = det(Nt)~ 1/2 MZt) ■ 

By Jacobi’s determinant formula dt det(N t ) = det ( A T t ) tr (dtN t Nf 1 ) we have 

dtdet(Nt)- 1 / 2 = -Jdet(Nty^tiidtNtNf 1 ) . 

We now use the Hamiltonian systems 

d t St = 9-HtSt , d t s; = SfH t 9 T , 

to differentiate the normalisation property b{StZQNt)*9l{StZt)Nt) = Id n . We obtain 


0 = d t N*N t ~* + iJV t *( S t Z 0 y (Ht - Ht)(S t Z 0 )Nt + N^ l d t N t 
= dtNtNf 1 - Z* Im H t Z t + N- ] d t N t , 


and by Proposition |2.5| 

tr (dtNtN^ 1 ) = 2 tr(Z f * Im H t Zt) = \ tr (Im H t {G~[ l — ill)) = \ tr(Im H t Gf 1 ) . 
It remains to write det(W) -1 ^ 2 = : an d to observe that 

dtPt = - j tr (Gf l Im H t ) , /3 0 = 0 . 


□ 
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The imaginary part Im Ht together with the symplectic metric Gt E Sp(n, R) determine 
the real-valued scalar 

A = - \ [ tr(G“ 1 1m H T )dr 

Jo 

and therefore the norm of the propagated coherent state, 

||C/(%o(^o)|| = ||e^Vo(A)||=e^ . 

The evolution of the symplectic metric Gt and the corresponding complex structure Jt = 
—QGt are governed by the following Riccati equations. 

Theorem 4.3 (Riccati equations). Let Lq and Lt = StLo be positive Lagrangian subspaces. 
Denote by Gt,Jt E Sp(n, R) the symplectic metric and the complex structure of Lt, respec¬ 
tively. Then, 


G t = Re H t QG t - G t D Re H t - Im H t - G t Q Im H t QG , 
j t = Q Re H t .J t - J t D R eH t + n Im H t + ./,< > Im H t J t . 

Proof. The equations of motion for Gt and Jt have been derived in |GS12( Theorem 3.3] using 
the Siegel half space and rational relations. The appendix [B] provides an alternative proof 
based on Lagrangian frames. □ 


4.4 Excited state propagation 

Next let us consider the propagation of first order excited states A**{1)<pq{Zq) for l E Lq. By 
Lemma |4.1| and Proposition |4.2| we obtain 

U(t)A'(l)<p 0 (Z 0 ) = U{t)A\l)U-\t)U{t)ip 0 {Z Q )=^A\S t l)MZt) ■ 


If St is a complex matrix, then Stl ^ L t = S)Lo so that A'(Stl) is not a creation operator 
associated with L t . We therefore use Proposition T3 and decompose 

Stl = TTL t Stl + TTZ Stl = TTLt.Stl + TJjStl , 


which leads to 


A*(S t l) = A\ir Lt S t l) + A\wr t S t l) = A\ir Lt S t l) - A(n Lt S t I ) . (4.6) 


Therefore, 

A\S t l)MZt) = A^(n Lt S t l)MZt) , 

since A(TTL t Stl)(po(Zt) = 0. The following Lemma and Theorem extend this line of argument 
to vectors of excited states. 

Lemma 4.4 (Ladder decomposition). Let L$ and Lt = StL q be positive Lagrangian subspaces 
and Zt E F n (L t ). Then, 

A\S t Z 0 ) = C* t A\Z t ) - DjA(Z t ) , 

where Ct = ^ZfD T StZ^ and D t = fZfLl T StZo are the unique matrices in C nxn so that 

S t Z 0 = ZtCt + Z t D t . 
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Proof. We apply the decomposition ( |4.6[ ) to the column vectors l\, ..., l n of Zq and obtain 

A^S t Z 0 ) = A\ir Lt S t Z 0 ) - A(n Lt S t Z 0 ) . 


Now we want to find Ct and Dt such that 

7 T Lt S t Z 0 = Z t Ct and 7r Lt S t Z 0 = Z t D t (4.7) 

because then (nL t StZo) = CfA^(Z t ) and A(iTL t StZo) = DjA(Z t ), which will give the result. 
We just multiply the equations in (4.7) from the left by ^ZfD T and use the normalisation 
of Zt, which gives 

Ct = l -Z* t n T ^ Lt S t Z 0 and D t = l -Z;Q T 7T Lt S t Z 0 . 

Now it remains to compute ^ZfLl T TT^ t = ^ZfLl T fZt.Z*Ll T = ^Z*Ll T , and this leads to the 
claimed expressions for Ct. and Dt. Adding the defining equations in (4.7), we finally obtain 

ZtCt + ZtDt = TTL t StZo + T^L t St Zq = StZo . 


□ 


Since A(Z t )(po(Z t ) = 0, the Lemma implies 


U(t)A\Z 0 )ip 0 (Zo) = e^A\S t Z 0 )MZt) = e A CjA*(Z t )<p 0 (Z t ) . 


(4.8) 


Let us next consider more highly excited states and expand 

C/(t)A^(Z 0 )(/? 0 (^o) with |t| > 1 

in terms of the orthonormal basis <p a (Zt ), ct G INq. Now the annihilation part of the decom¬ 
position 

A\S t Z 0 ) = C* t A\Zt) - DjA(Z t ) 

becomes more visible and we encounter commutators between CfA^(Zt) and DjA(Zt). In 
this situation, the term handling is facilitated by multivariate polynomial recursions that are 
governed by a complex symmetric matrix. 

Theorem 4.5 (Excited state evolution). Let Lq and Lt = S)Lo be positive Lagrangian sub¬ 
spaces. Consider Zt G F n (Lt) so that Zt = StZ^Nt for a Hermitian positive definite matrix 
N t G C nxn , and denote by Gt G Sp(n, H) the symplectic metric of L t . Define 

M t = \(StZ 0 ) T Gt(StZ 0 ) 

and the polynomials q a {x), x G C n , a G INq, via the recursion relation 


q 0 (x) = 1, q a+ej (x) = xjq a (x) - ej ■ M t Vq a (x) , j = l,...,n. 


Then, we have for any a G INq 


U (t)ip a (Zo) 




q a (N t A\Z t ))MZt ) • 


(4.9) 


(4.10) 
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Before entering the proof, we briefly examine the special case of Hermitian time evolution. 
In this case St is real and we can choose Z t = St.Zo- Then, N t = Id n and f3t = 0, and 
Proposition 


2.5 


implies Mt = GtZt = \Z( QZt = 0. In summary, 


1 


U(t)<p a (Z 0 ) = -—A f a (Z t )ipo(Z t ) = <p a (Z t ) 


t) , 


which is of course also directly implied by U{t)A\Zo)TJ(t ) = A^(StZo), see Lemma 4.1 or 
|Hag80| and jH ag98| . In the more general non-Hermitian case we observe that the time 
evolution activates lower order states. Equation (4.10) can be interpreted as an expansion of 
the propagated state into the basis defined by Z t , 

U(t)<p a (Z 0 ) = e* E ^-kNkX^t) i 

\k\<\a\ 


where the coefficients a & can be computed in terms of Nt and the polynomial q a . 

It is worth emphasising the prominent role played by the matrices Nt and Mt. All the 
information about the effects of the non-Hermiticity on the propagation are encoded in those 
two matrices. 


Proof. We have by Lemma 



ejA\S t Z 0 ) = eJC* t A\Z t ) - ejD?A(Z t ) =: Uj - Vj . 


Then, 

n 

Al(S t Z 0 ) = H(u j -v j )N , 

3 = 1 

and in particular A ] a+e . (S t Z 0 ) = ( Uj - Vj)Ai(S t Z 0 ). If we apply At(S t Z 0 ) to ipo{Z t ) <E I{L t ), 
then we can use that Vjipo(Zt) = 0. We then commute all Vj to the right of the Uj and obtain 
that 

A' a (S t Zo)<po(Z t ) = q a (ui,--- , u n )<po(Z t ) , 

where q a {x ) is a polynomial in n variables. Our aim is now to derive a recursion relation 
for q a {x). Let us define a matrix M = M t G C nxn by 

Mij := [vi,Uj\ . 

Then we have VjvJf = [dj. u^] +v%Vj = Mj^ku^ 1 +u^Vj and for any polynomial p(u\, • • • , u n ) 

[vj,p(ui, ■ ■ ■ ,u n )\ = (ejM t \7p)(u i, ■■■ ,u n ) . 

We therefore find 

Vjq a (ui, ■■■ ,u n ) = ejMtVq a (ui, ■■■ ,u n ) + q a (ui , ■ ■ ■ ,u n )vj . 

Using q a+ej (ui, • • • , u n ) = (uj ~ Vj)q a (ui, ■■■ ,u n ) and Vj(p 0 (Zt) = 0 we get 

q a +ej(ui, ■ ■ ■ ,u n )ipo(Z t ) = Ujq a (ui, ■ ■ ■ ,u n )<po(Z t )-ejMtVq a (ui,--- ,u n )po(Z t ) , 
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which is the recursion relation ( |4.9[ ). It remains to compute Mf. We obtain from Lemma 3.2 
that 

[vi,Uj] = [A(Z t D t ei), (ZtCtej)} = \{Z t D t ef) T = eJ{DjCt) ej . 


Hence, by Lemma |4.4| and Proposition |2.5| 


M t = DfCt = IZ^S^nZtZf^StZo = l(S t Z 0 ) T (Gt + ifi)(S t Zo) 

= \(S t Z 0 fG t (S t Z 0 ) . 

We observe that M t = Mj and notice that (tti, • • • ,u n ) T = CfA\Z t ). Moreover, 

Ct = ^Z*(S^n T S t )Z 0 = l -N t z*n T z Q = N t , 

since St is symplectic and Zq normalised. □ 

Applying a linear combination of powers of A'(Zt) to the normalised Gaussian ipo(Zt;x) 
produces multivariate polynomials in x that can be described by a recursion relation of the 
type encountered above. 

Corollary 4.6 (Polynomial prefactor). Let Lq and Lt := S)Lo be positive Lagrangian sub¬ 
spaces. Let Zt E F n (Lt) so that Zt = StZoNt for a Hermitian positive definite matrix 
N t E C nxn and set 



Denote by Gt E Sp(n, K) the symplectic metric of L t . Define 

M t = \(StZ 0 ) T Gt(S t Z 0 ) and M t = M t + N t Qf l Q t N t ■ 

We then have for any a E INq 

U(t)(p a (Z 0 ]x) =-^ PafyJlNtQ^x) <p 0 (Z t ; x) (4.11) 

where the polynomials p a (x), x E C n , satisfy the recursion relation 

Po(x) = 1 , p a+ej {x) = XjPa{x) - ej ■ M t Vp a (x) , j = 1,... ,n . 

Proof. We first compute 

N t A\Z t )<p 0 (Z t -x) = ^NtiPfQt-QlPjQf'xMZux) 

= \fl N tQt lx( Po(Z t ;x) 

= yhpo{z t -,x) , 

with yt = Nt.Q~f 1 x and where we have used the normalisation ZfLlZt = Q\Pt~ PfQt = 
2ild n . This motivates the ansatz 

q a (N t A ] (Z t ))ip 0 (Z t ; x) =: p a (y t )po{Z t \x) . 
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The gradient formula of Lemma |C.2| implies 

Pa+ ej (yt)vo{Z t -,x) = ej ■ N t A ] (Zt) (p a (y t )ip 0 (Z t -x)) - ej ■ M t (a j p a - e .(yt)) , . =1 (po(Zt;x) . 
We compute 

N t A^(Zt)(p a (yt)<Po(Zt;x)) = p a (yt)NtA\Z t )ipo(Z t -,x) - {Z t \x)N t Q* t pp a {yt) 

= ytp a {yt)vo{z t ; x) - N t Q* t Qi T Nj(yp Q )(y t )(p 0 (z t -,x) 

so that 

Pa+ ej (yt) = ytjPaiyt ) - ej ■ (M t + N t Q* t Q~ T N^)(a j p a ^ ej (y t )) 1 - =1 . 

Since QtQt is real symmetric, we have Q* t Qf T = Qf l Qt and M t is symmetric. □ 


4.5 Dynamics of the center 

Repeating the calculations of Lemma |4. 1| , the dynamics of the centered ladder operators read 

U(t)A{l,z)U-\t) = A(S t l, S t z) , (4.12) 

U{t)A\l,z)Xr\t) = A\S t l,S t z) (4.13) 


for all l, z E C n © C n . We assume that Lq and L t = StL$ are positive Lagrangian subspaces 
and consider the complex structure Jt E Sp(n, R) of the Lagrangian L t . We then know by 
Theorem 3.12 that a real projection of the center Stz does not change the ladder operator if 
we parametrise by the Lagrangian Lt , that is, 


A(S t l,S t z) = A(S t l,P Jt (S t z)) 


for all l E Lq and z E C n © C n . The dynamics of the projected center are easily inferred from 
the Riccati equations for the complex structure ,Jt- They reflect the metriplectic structure of 
equation (4.4) on the finite dimensional level. 


Corollary 4.7 (Projected dynamics). Let Lq and Lt = StLo be positive Lagrangian sub¬ 
spaces. Denote by Gt, Jt £ Sp(n, R) the symplectic metric and the complex structure of Lt, 
respectively. Let zq E R n © R". Then, zt '■= Pj t (StZQ ) E El™ © R n satisfies 


z t = D Re HfZt + G t 1 Im H t zt . 


(4.14) 


Proof. We differentiate zt 


R e(Stz) + Jt Im(S t z) so that Theorem 4.3 implies 


z t = Re(f 1H t S t z) + j t Im(St^) + J t Im(S lH t S t z) 

= Re Jit R e(Stz) — Im Ht Im(S't^) 

+ (n Re H t J t - J t D ReHt + Q Im H t + J t D Im H t Jt ) Im(S't^) 
+ J t Ll Im H t Re(Stz) + JtLl R e(H t ) Im(Stz) 

= O Re HtZt + Jfil Im /// zt . 


Moreover, LlJ t = Gt gives Jfi-l = ii'GfiJ- = G t 1 . 


□ 


The time evolution of coherent states with real projected center resembles the one of 
Hermitian dynamics, however, with a phase factor determined by the action integral of the 
Hamiltonian Tit along the real projected trajectory. 
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Proposition 4.8 (Coherent state evolution). Let Lq C C n © C n be a positive Lagrangian 
subspace, Zq E F d (Lq) and zq E Et n ©R”. Let the coherent state pq(Zq,zq) be given by (3.8). 
If the Lagrangian Lt = SjLo is positive for t E [0, T[, then 

U(t)MZo,z 0 ) = <p 0 (S t Z 0 ,z t ) = ei^H^o (Z t ,z t ) 

for all t E [0,T[, where zt =: ( pt,qt ) £ is defined by (4.14), fit is the factor derived 

in Proposition \4-I\ and 

a t (z 0 ):= (q T ■ p T - U T (z T )) dr (4.15) 

Jo 

denotes the associated action integral of the Hamiltonian Ht along zt ■ 

Proof. We use the time evolved lowering operator ( |4.12[ ) with the real center zt and obtain 

U(t)A(Z 0 ,z 0 )<po(Z 0 ,z 0 ) = U(t)A(Z 0 ,z 0 )U~ 1 (t)U(t)po(Z 0 ,zo) 

= A(S t Z 0 , z t )U(t)ip 0 (Z 0 , z 0 ) . 

Therefore, ipo(Zo,zo) £ I(Lq,zq) implies U(t)(po(Zo, zo) £ I(L t ,Zt), and hence there exists 
ct E C with U(t)ipo(Zo, zo) = ctp>o(St.Zo, zfi) =: c±tp o(t). It remains to determine q. We 
denote 

'Pt 

K Qt, 


H t = 


Hpp Hpq 
Hqp Hqq 


S t z 0 = 

Computing iedt(ctpo(t)) we obtain 

iect/ct + ie (d t det(Q t )“ 1/2 j det(Qi) 1/2 + \ed t (^(x - q t ) ■ B t (x - q t ) + fpt ■ (x - q t )) 

times ctpo(t). We sort this second order polynomial in powers of ( x—qt ) and keep the constant 
terms, that is, 

iect/ct - f tr (d t Q t Qf l ) + Pt ■ <lt , (4.16) 

where we have used Jacobi’s determinant formula (dt det Qt)/ det Qt = tr (dtQtQfi 1 )- Next 
we compute 

BtPo(t) = \(z • H t z)ip 0 (t) = ■ ((HppB t (x - q t ) + Hp P p t + H pq x) ip 0 (t)) 

© 2 ® ‘ (HqpBfix qt) + HqpPt © HqqX)) ipo(t) 

Therefore HtCtPoit) is a second order polynomial in (x — qt) times ctp>o(t), and the constant 
terms amount to 

| tr (HppPtQfi 1 + H pq ) + Ht(z t ) . (4.17) 


Since dtQt = H pq Q t + Hp p P t and 8tQtQ t = H pq + H pp P t Q t , the matching of the terms in 
(4.16) and (4.17) gives 

i ect/ctPpt ■ qt = Bt{z t ) , 

which is solved by the exponential of the action integtral q = □ 


Our previous results on excited state propagation, that is, 
describe the time evolution of 


Theorem 4.5 and Corollary |4.6[ 


U(t)p a (Z 0 ,z 0 ) , a E INq" - , 

for the case zq = 0 in terms of multivariate polynomials. Essentially, these results stay the 
same when considering nonzero zq E IR^ © We only have to record the evolution of the 
center and add the corresponding action integral. 
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Theorem 4.9 (Excited state evolution). Let Lq and Lt 
spaces. Let zo £ H n © R n and Z t E F n (L t ) so that Z t 
definite matrix Nt E C" xn . Set 


Z t = 



StLo be positive Lagrangian sub - 
StZoNt for a Hermitian positive 


and denote by Gt E Sp(n,R) the symplectic metric of L t . Define 


M t = \{S t Z 0 ) T G t (S t Z 0 ) and M t = M t + N t Qfi l Q t N t . 


Then, we have for any a E INq 


U{t)(p a (Z 0 ,z 0 ;x) = 


e^ 


^ at (z 0 )+Pt 


q a (N t A\Z t , z t )) tpo(Z t , zp, x) 


e£ 


^a t (z 0 )+Pt 


Pa 


i N t Q t (x-q t )) To{Z t , z t ; x) 


where zt = ( pt,qt) £ R n © R n is defined by ( |4. 14 ) and afizfi) is the action integral \f-15 ) of 
Tit. along the trajectory zt■ The polynomials q a {x) = r a (x-,M t ) and p Q (x) = r a (x; M t ) satisfy 
the recursion relations 


r 0 (x;M) = 1 , 


r a+ej (x: M) = Xjr a (x; M) - ej ■ MVr a (x; M) , j = 1,..., n , 


with M = Mt and M = Alt, respectively. 


The time evolution of almost all the constitutive elements of Theorem [4JJI can be described 


by ordinary differential equations: First, there is the Riccati equation of Theorem 4.3 for the 
symplectic metric Gt, that can be solved together with the equation for the loss or gain 
parameter fit, 


dtfit = -3 tr (G t 1 Im H t ) , ,% = 0 . 

Second, there is the metricplectic equation of Corollary |4.7| for the real center zt, together 
with the corresponding action integral at(zo). Finally, for the normalised Lagrangian frame 
Z t = St.ZoNt, we find the equation 


dtZ t = QHtZ t + Z t Nfi x d t N t . 


which contains the time derivative of the normalising matrix Nt. We will illustrate in the 
following section how one can determine Nt for an explicit example. 


5 The Davies—Swanson oscillator 

As an example we investigate the dynamics of a one-dimensional quadratic non-Hermitian 
Hamiltonian, the Davies-Swanson oscillator 

= ^(j5 2 + q 2 ) - ^{pq + qp) = \ Op[z • Hz] 

defined by the complex symmetric matrix 
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whose imaginary part is a real symmetric matrix with eigenvalues ±5. For this particular 
Hamiltonian the spectrum and transition elements have been computed jDav99at ISwa04j as 
well as the dynamics of coherent states [OKR.S14] . It is our aim here to complement the 
picture by propagating excited wavepackets. 


5.1 One-dimensional systems 


For one-dimensional systems, the results of Theorem 4.5 simplify, since the normalisation 
of Lagrangian frames just involves the inversion of a positive real number. Starting with a 
positive Lagrangian subspace Lq = spanjZo} spanned by a normalised vector Iq E C © C, we 
set 

^ j with nf 2 = h(Stlo, StJo) > 0 


k ■— Stlont — 


to obtain a normalised Lagrangian frame It E C © C of the time evolved subspace Lt = StL q. 
In order to describe the propagation of excited wavepackets we use 

m t = dtct = n 2 h(S t lo , S t l o) . 

For notational convenience, we restrict ourselves to the case zq = 0. For non-vanishing 
centers zq E R © R, there is an additional multiplicative factor due to the complex-valued 
action integral ott(zo). According to Proposition |4.2| and equation (4.8), we obtain 

U(t)ip 0 {l 0 ) = e^ipoih) , 

U(t)ipi(lo) = e^ntipiUt) , 

U{t)tf2(lo) = ^{ntWik) ~ mnpoilt)) 

for the coherent and the first two excited state, respectively. For the whole orthonormal basis 
Theorem 4.5 provides 




U(t)ipk(lo ) = -^qk{n t A ] {l t ))yo{lt 


k € 1N 0 , 


where the univariate polynomials q & satisfy the recursion relation 

go 0 ) = 1 , q k+ i(x) = xq k (x) - m t q' k (x) , t£l 0 


(5.1) 


(5.2) 


These polynomials are Hermite polynomials with time dependent scaling according to the 
complex number mt. Using the monomial expansion of these polynomials, we can rewrite the 
expansion in 0 explicitly in terms of the propagated basis functions (pk{h )■ 

Corollary 5.1 (Explicit expansion). Let Iq G C 0 C so that L$ = spanj/o} and L t = <S)Lo 
are positive Lagrangian subspaces. Let nt > 0 so that It = Stlont is normalised according to 
h(lt,lt) — 1. Then, for all k E INq, 


U(t)tp k (l 0 ) = ^2 

3=0 


afajUf <p>j ( It 


where a k j E C are the coefficients of the monomial expansion q k (x) 
polynomials q k (x) defined by the recursion relation (5.2). 


Y]j =o a kjX j of the 
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Figure 1: Time evolution of \\U(t)ipk(lo)\\ according to Corollary 5.1 with Iq = (1, —i). The 
color coding is blue for k = 0, red for k = 1 and yellow for k = 2. The parameters of the 
oscillator are chosen as cuo = 1 and 5 = 0.5. In all three cases the norm considerably depart 
from unity, the larger k is, the stronger they deviate. 


Proof. Since for all j, we have 

k 


U(t)ip k (l 0 ) = -j=J2 a kjrit A ](lt)Po(lt) = e0t '52y jj a kjn J t (fj(lt) 
vfc- j= o j=o 


□ 


5.2 Norm evolution for the Davies—Swanson oscillator 

Applying the one-dinrensional formulas to our particular example, the Davies-Swanson oscil¬ 
lator, we start by examining the classical Hamiltonian system 

S t = VHS t , S 0 = Id 2 ■ 

Its solution St = exp(tQH) exists for all times t € It. Setting u 2 := Uq + 5 2 , we observe 
(QH) 2 = —cu 2 Id 2 and consequently 

(PlH) 2k = (-l)V fc Id 2 , ( QH) 2k+1 = (-1 ) k uj 2k QH , k > 0 . 

Therefore, 

00 + 2k 00 j.2fc+l 

k — 0 k — 0 

= cos(tu;)Id 2 + ^ sin (tufjVLH . 
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This formula for St allows to explicity compute the time-intervals for which a particular initial 
Lagrangian subspace stays positive. 


Lemma 5.2 (Positive Lagrangian subspace). Let Lq = span{Zo} with Iq = (1, —i) 6 C © C 
and consider L t = StL®. If ui o > <5, then L t is a positive Lagrangian subspace for all i £ 1. 
Otherwise, L t is positive for t £ [0, T[ with 


T := — arccos 
2 c j 



Proof. We first compute 

SfLiSt = (cos(ftu)Id2 + ^ sin(tuj)LlHy (cos(tu;)Id2 + ^ sm(tu)QH) 
= cos 2 (tui)Ll + J 2 sin 2 (tui)HflH — 2 \ cos (tu) sin(ta;) Irn H . 


and 


HUH = 


2\5vjq 
u 2 - 6 2 


S 2 - ^ 

— 2\6ojq ) ’ 


Irn H = 


This implies for all normalised vectors / 6 C © C with hfl, l ) 


0 -A 
-5 0 ) ' 

2 l*Ll T l = 1 that 


h(Stl,Sfl) = cos 2 (tui ) — 2^2 sin 2 (tui)l* HQ HI — I cos(tuj) sin(tuj)l* Im HI . 
With l = (h,l 2 ) 

Pimm = 2i5coo(\h\ 2 - \l 2 \ 2 ) - 2 i(wg - 5 2 ) Im(Ii/ 2 ) , 

Z*ImH7 = — 2J Re(7iZ 2 ) , 


In one dimensional systems the normalisation of l is equivalent to Im(ZiZ 2 ) = 1, so we can 
replace imaginary part in the equation above. However, there is no relation between l\ and 
l 2 in general, so we cannot simplify this further. 

For the particular vector Iq = (1, —i) we obtain 


Wn - 5 2 


5 2 


h(S t lo, S t l o) = cos (to;) + ^ siir(tw) = 1 -, (1 - cos(2 too)) . 

uj o + cP 

This function is positive for all t 6 R, if uiq > 6. Otherwise positivity holds on [0,T[. 


□ 


We work for times t so that the Lagrangian subspace L t = 5)Lo is positive and consider 
the normalisation factor m > 0 with 


n 


-2 


s 2 

= h(S t lo, Stl o) = 1-^ (1 - cos(2tui )), l 0 = (1, -i) . 

oi- 


Then, 


It Stlont — nt 


cos(fcj) + i sin (tuj) 
— icos(tcj) + sin(icT 


€ F n (L t ) 


is a normalised Lagrangian frame and the symplectic metric of Lt reads 
Gt = Q t Re(ltlt)Ll = n 2 ld 2 + 2 ^n 2 sin(Zo;) 


— ^sin (tui) cos(tui) 

cos(tui ) ^ sin (tu) 
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An elementary calculation yields 


2 A 2 

ti^G^ImiL) =-n 2 sin(2Lj) 

UJ 


so that 


and 


Pt = J tr(G T 1 ImH)dr = - \ ^ln^^-^ - Inf ^ + cos(2 tot) 


?P* = u 1//2 (wq + <5 2 cos(2o;t)) 


1/4 


For the polynomial recursion (5.2) we also have to compute 

m t = n\ h(S t l 0 ,Stk) . 


Repeating a part of the calculations of the proof of Lemma 5.2 we obtain for alH E C © C 


so that 


FHQHl = 2\5 uq (Zi - l 2 ) and f Im HI = -25 hl 2 
25 


ZO 2 . / /w 0 . / \ . , ,\ 

mt = — n t sm(ru;) I — sm(fu;J + icosfiu;) I . 

Having derived explicit formulas for the time evolution of the parameters, we now use Corol¬ 
lary 5.1 and compute the norm of the coherent state and the first two excited states as 

\\U(t)(p 0 (l 0 )\\ = e A , \\U(t)ipi(l 0 )\\ = e Pt n t , \\U(t)<p 2 (lo)\\ = + \\m t \ 2 . 

As expected, all three norms considerably depart from unity, the more highly excited the 


state, the stronger the deviation, see Figure 5.1 


A Weyl calculus 

Let us recall a few standard results about products and Weyl quantisation, see [ CR12 , Chap¬ 
ter 2] for background. We consider smooth phase space functions a, b so that 


Op[a]V’0) = (2vre) n [ a(£, \{x + y))e^ x v ' > ip(y)d^dy 

Jr 2 " 

together with the compositions Op [a] Op [6] and Op [b] Op [a] are well-defined linear operators 
on dense subsets of L 2 (R n ). The symbol of the operator product is the so-called Moyal 
product of a and b, 

Op [a] Op [b] = Op[aj|6] . 

If one of the two symbols a or b is a polynomial of degree < 2, then 


ajjft = ab + — Va • HV6 — — tr ( D 2 aQD 2 bQ T ) , 
2 8 


where V = V Pj(? and 



—Id A 

0 J 


G R 


2nx2n 
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Consequently, the commutator can be written as 


[Op[a], Op[6]] = ie Op[Va • nV6] . 

In particular, the canonical commutation relations can be quickly verified as 

O' 


(A.l) 


\qj,p k \ = i eVqj • Wp k = ie 


e-k 


= ieSj 


1 jk ■ 


Another application of the product rule yields that the Weyl quantisation of a symmetric 
quadratic form equals the quadratic form in z. 

Lemma A.l (Quadratic symbol). We consider 

H=(H pP H™) € € 2nx2n ' 

\-Hqp tiqq) 

Then, Op[z • Hz] = z • Hz + tr (H qp — H pq ). In particular, 

Op [2; ■ Hz] = z ■ Hz , if H = H T . 

Proof. We compute 


2 n 


■ Hz = Op [z ■ Hz] + j ^2 Vzj • QH jk Vzk = Op[z • Hz] + ^ tr (-Hpq + H qp ) , 


j,k =1 


since 1 ■ Hcfc) — H\ n • • • H n 2 n T H n \ T • • • T H 2 n,n — tr( H pq T H qp ) 


□ 


B Dynamics of the metric and the complex structure 


We provide a Lagrangian frame’s proof for Theorem |4.3[ that states the Riccati equations for 
the symplectic metric and the complex structure of the positive Lagrangian Lt = StL 0, that 
is, 

G t = Re H t QG t - G t Q Re H t - Im H t - G t n Im H t QG , 

.// = ( > Re It / •// — .// Q Re /// + II Im It / + Jfl Im /// ,7/ . 


Proof. We only work for Jt , since Gt = 
matrix jVj € C nxn so that Zt = StZoNt 
the proof of Propositon 4.2 we obtain 


fUt. Let Zo e F n (Lo ) and consider an invertible 
6 F n (L t ). We then have Jt = — R e(Z t Zf)Tl. As in 


0 = 5tIV*IV t -* + ilVt(5tZo)*(//t-^t)(5tZo)At + A t - 1 5tIVt 
= dtIV; JV t "* - Z* t Im H t Zt + N^dtNt . 


Next we differentiate Zt so that 

Z t = nH t S t Z 0 N t + S t Z 0 d t N t = QH t Z t + Z t N^ l d t N t . 

Therefore, 

5t(ZtZ*) = QHtZtZ? + ZtN^dtNtZ? + ZtZ;H t n T + Z t d t N t *N t ~*Z* 

= nH t z t z* + z t z;H t n T + z t z t * im LTtZtz;. 
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Since Im (ZfZ*) = —H, we then have 

d, R e{Z t Z*) = <> Re //, R e{Z t Z*) + O Im //,<> - Re(ZiZ*) Re //,Q 
+ Re(Z 4 Z;) Im Re(Z t Z t *) 

and the claimed equation .// = D Re /// ,// + 11 Im f// — J/Q Re f// + .//11 Im H/ J/ . □ 


C Multivariate polynomials 


Analysing Hagedorn wave packets and their dynamics, we have encountered multivariate 
polynomials generated by the following type of recursion relation. 

Definition C.l (Polynomial recursion). Let M £ C nxn be symmetric and c £ C. We define 
a set of multivariate polynomials p a {x ) by the recursion relation 


Po(x ) = c , 

with x £ C n and a £ INq . 


(p a+ej .(x))" 1 = xp a (x ) - MVp a (x) 


(C.l) 


Together with c = 1 the matrix M = 0 generates the monomials p a (x) = x a , while the 
identity matrix M = Id determines tensor products of simple Hermite polynomials. 

If Q £ C nxn is the lower block of 

' P \ 


Z = 


Q 


then the unitary matrix M = Q l Q generates the polynomial prefactor of the Hagedorn wave 
packets, that is, 


1 


p a (z-,x) = y=pa (yi yj <Po(Z-,x) with y = Q x. 

All the polynomials sequences of Definition |C.l are multivariate versions of the so-called 
Appell sequences. In the univariate setting, a polynomial sequence p n (x), n £ INo, is called an 
Appell sequence, if p' n (x) = np n -\{x) for all n > 1. The Hermite polynomials are prominent 
examples. In several dimensions this property generalises to the following gradient formula, 
which is due to [ DKT15 J. 

Lemma C.2 (Gradient formula). Let M £ C nxn be symmetric and c £ C. The polynomials 
defined by the recursion relation (C.l) satisfy 


V x.Pa(x ) = (ajPa-e^x))™ 


for all a £ INq and x £ C n . 


Proof. We argue by induction and assume that the gradient formula holds for a fixed a £ INq . 
Differentiating the recursion relation, we get 

dkPa+ej = 5 k jPa + Xjd k p a ~ ej ■ MV d k p a = 5 k jPa + Ot k (xjp a -e k ~ ej ■ MX7p a -e k ) 

— 3kjPa T &kPa+ej—e k — (® T S-fikPa+ej—e k • 


□ 
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